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Earth rotation prevents exact solid body rotation of fluids in the 
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Abstract -We report direct evidence of a secondary flow excited by the Earth rotation in a water- 
filled spherical container spinning at constant rotation rate. This so-called tilt- over flow essentially 
consists in a rotation around an axis which is slightly tilted with respect to the rotation axis of 
the sphere. In the astrophysical context, it corresponds to the flow in the liquid cores of planets 
forced by precession of the planet rotation axis, and it has been proposed to contribute to the 
generation of planetary magnetic fields. We detect this weak secondary flow using a particle image 
velocimetry system mounted in the rotating frame. This secondary flow consists in a weak rotation, 
thousand times smaller than the sphere rotation, around a horizontal axis which is stationary in 
the laboratory frame. Its amplitude and orientation are in quantitative agreement with the theory 
of the tilt-over flow excited by precession. These results show that setting a fluid in a perfect solid 
body rotation in a laboratory experiment is impossible — unless tilting the rotation axis of the 
experiment parallel to the Earth rotation axis. 
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Introduction. — There are few examples of fluid me- 
chanics experiments at the laboratory scale in which the 
Earth's Coriolis force has a measurable influence. Such 
experiments may be considered as fluid analogues to the 
Foucault pendulum. The most popular instance is cer- 
tainly the drain of a bathtube vortex pQ . Although this is 
the subject of common misconception, it is actually pos- 
sible to detect the influence of the Earth's rotation on 
the vortex, but only under extremely careful experimental 
conditions, far from the everyday experience [2]. Ther- 
mal convection is another example, in which a slow drift 
qf the large-scale flow due to the Earth rotation has been 
detected in very controlled systems [3JH]. 

In this letter we describe an experiment which may be 
considered as the most simple fluid Foucault pendulum: 
it consists in a volume of water enclosed in a spherical 
container spinning at constant rotation rate (fig- [I])- 
After a transient known as spin- up, the water is expected 
to rotate as a solid body at the same rate Oq [5]. The 
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timescale for this spin-up is classically given by the Ek- 
man time te = R (V^o) -1 ^ 2 , where R is the radius of the 
sphere and v the kinematic viscosity of the fluid. For a typ- 
ical laboratory experiment using water, this timescale is 
usually of order of a minute, so after a few tens of minutes 
a perfect solid-body rotation should be reached, with the 
fluid exactly at rest in the frame of the container. If this 
simple experiment is performed on Earth, it is expected 
that the Earth rotation could prevent from reaching this 
idealized solid rotation state [6j[7] . A weak secondary flow, 
known as tilt- over flow [318 ,9 , is induced by the preces- 
sion of the rotation vector of the container by the Earth 
rotation vector Cl p . Seen from the laboratory frame of ref- 
erence, the fluid particles rotating at velocity uo = Flo x r 
experience a Coriolis force per unit mass f c = —2Cl p x uo- 
This Coriolis force disturbs the fluid particles periodically 
at frequency Qq, and tends to deflect their trajectory to- 
wards the plane normal to fl p . 

Precession driven flows in spherical or spheroidal con- 
tainers and in spheroidal shells have received considerable 
interest since Poincare [10 , because of their importance to 
geophysical and astrophysical flows [8,9j. In the case of the 
Earth, rotating with a period Tq ~ 1 day, the precession 
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Fig. 1: (Color online) Sketch of the rotating platform and the 
water filled sphere. (ei,e2,es) is a Cartesian coordinates sys- 
tem attached to the laboratory frame. The platform is rotating 
at flo = Qo ©3 in the laboratory. Q p is the Earth rotation vec- 
tor at the latitude A = 48.70° of the laboratory, oj is the 
rotation vector of the tilt-over flow in the bulk. The rotation 
vectors are not to scale. 



of its rotation axis, at a period T p ~ 26 000 years, could 
produce large excursions of the rotation axis of the liquid 
core [11] . Precession driven flows have also been proposed 
by Malkus [9] to contribute to the generation of planetary 
magnetic fields, which has been later confirmed by Ker- 
swell [12] and Tilgner [13]. Kida [14] recently proposed a 
complete solution for the flow in a rapidly rotating sphere 
under weak precession, including a detailed analysis of the 
conical shear layers detached from the critical latitudes. 

First evidence of a tilt-over flow excited by the Earth 
rotation in a laboratory experiment has been reported 
by Vanyo and Dunn [6], using visualizations by dyes 
and buoyant tracers, but without quantitative determi- 
nation of the tilt-over flow properties. Recently, Triana 
et al [7] obtained indirect evidence of this effect, from 
one-dimensional velocity profiles in a rotating water-filled 
spherical shell, 3 m in diameter, containing an inner co- 
rotating sphere. However, no quantitative agreement with 
the theory of Busse [8J could be obtained in their experi- 
ment. 

Based on the same idea, we provide in this letter, by 
means of particle image velocimetry measurements (PIV), 
the first direct visualization of the precession flow driven 
by the Earth rotation in a sphere rotating in the labora- 
tory. These measurements are a technical challenge, be- 
cause of weakness of the velocity signal of this tilt-over 
flow (the fluid rotation axis is tilted by less than 0.2° with 
respect to the sphere rotation axis). A quantitative agree- 
ment with the theory of Busse is demonstrated, both for 
the magnitude and the orientation of the secondary circu- 
lation. 



Physical origin of the tilt-over flow. — 

Poincare JTU] first analyzed the precession flow in a 
sphere in the singular case of a perfect fluid. He showed 
that the inviscid solution consists in a solid-body rotation 
around an axis parallel to f2 p , but of undefined amplitude. 
In the presence of weak viscosity, far from the boundaries, 
the tilt-over flow may still be described as a solid-body 
rotation, with a rotation vector uj tilted with respect to 
ffco? and stationary in the precessing frame (the laboratory 
frame here). We note in the following uj' = up — Qq the 
rotation vector of the fluid in the bulk measured in the 
rotating frame. 

Remarkably, the presence of viscosity, even weak, dras- 
tically changes the rotation vector of the fluid u compared 
to the inviscid solution of Poincare. The orientation and 
amplitude of cv for a viscous fluid are now non trivial func- 
tions of the Poincare number Q p /Qq and of the Ekman 
number E = v/(Q R 2 ). In the limit Q p /Q < VE <C 1, 
the rotation vector uj is almost equal to fio, and the small 
correction uj r is almost normal to CIq. This tilt-over flow 
has been described by Busse [8 as one among a dense 
family of inertial modes, of eigenfrequency given by Qo 
(see Ref. [5] for a general description of inertial modes in 
a sphere). When forced by precession, the magnitude uo' 
of this tilt-over flow can be determined by a simple bal- 
ance between the Coriolis torque (in the bulk) and the 
viscous torque (at the surface of the container). The Cori- 
olis torque is of order Y c ~ pR 4 f c ~ pR 5 Q p Qo cos A, with 
p the fluid density and cos A = x ^ p \/^to^p- The 
viscous stress is given by a ~ pvAu/5, where Au ~ uo'R 
is the small velocity jump between the container wall and 
the fluid bulk, and 6 = (v/Qq) 1 ^ 2 is the thickness of the 
Ekman boundary layer. The resulting viscous torque is of 
order T v ~ R 3 a ~ pvuu'R^/S. Balancing the two torques 
gives the simple relation 



E~ 1/2 n p cos A. 



(1) 



Although very weak, this tilt-over correction may be sig- 
nificantly larger than the Earth rotation rate in a typical 
laboratory experiment where E <C 1. 

Experimental Setup. — The experimental setup, 
sketched in fig.[U consists in a spherical glass tank, of inner 
radius R = 115 ± 0.25 mm, filled with water and mounted 
on the center of a precision rotating turntable of 2 m in 
diameter. We use two Cartesian coordinate systems, both 
with origin at the center of the sphere: (i) (ei, e2, 63), at- 
tached to the laboratory reference frame (fig. [1]), with ei 
pointing to East, e2 pointing to North and e3 along the 
vertical; (ii) (e^, e y , e z ), attached to the rotating platform 
(fig. [2]), with e z = e3, in which the measurements are 
performed. 

The platform is rotating in the laboratory frame with 
a rotation vector Qq = ^o e 3- The angular velocity Qq is 
varied between 2 and 16 rpm, with temporal fluctuations 
less than ±5 x 10 -4 . The Ekman number E = v/(QoR 2 ) 
varies between 3.6 x 10 -4 and 4.6xl0 _5 in this range of Qq. 
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Fig. 2: (Color online) Schematic view of the cubic water 
tank containing the 115 mm radius glass sphere, mounted 
together on the rotating platform. PIV measurements are 
achieved in off-centered vertical and horizontal planes, located 
at ?/mes = +22 mm and z mea — +22 mm, using a corotating 
laser sheet and a camera aiming normally at it. The angle 
6 {t) between the images and the North direction is determined 
using a continuous laser beam aligned along the North- South 
orientation and crossing the rotation axis. 

In fig.HJ the rotation vector of the Earth fl p is also shown, 
for the latitude A = 48.70° of our laboratory in Orsay. The 
relative scale of the vectors Qo and fl p is obviously not 
realistic in this figure: the Earth rotation rate is Vt p ~ 
6.9 x 10 -4 rpm ~ 27r/(l day), which yields a Poincare 
number Q p /Qo ranging from 3.5 x 1(T 4 to 4.3 x 10~ 5 . 

After the start of the platform rotation, we wait at least 
r w = 2 hours before data acquisition in order to reach a 
stationary regime. This waiting time represents at least 
30 te, where te = R (z^o) -1 ^ 2 is the Ekman spin- up time. 
This indicates that the solid-body rotation state should 
be reached, apart from precession effects, with a relative 
precision better than exp(— t w /te) — 10 -13 . 

Velocity fields are measured in the rotating frame us- 
ing a two-dimensional PIV system [15] mounted on the 
rotating platform, in either a vertical (e x ,e z ) or a hori- 
zontal (e x ,e y ) plane (fig. [2]). These measurement planes 
are off-centered, at y mes /R — z mes /R 0.19 (see fig. [2]), 
in order to get better insight in the spatial structure of 
the flow. Optical distortions are reduced by immersing 
the glass sphere in a square glass tank of 300 mm side also 
filled with water. The distortion is found less than 5% for 
r < 0.9R. The fluid is seeded with 10 /im tracer particles, 
and illuminated by a corotating laser sheet generated by 
a 140 mJ Nd:YAG pulsed laser. For both horizontal and 
vertical measurements, the sphere cross-section is imaged 
with a high resolution 2048 x 2048 pixels camera aiming 
normally at the laser sheet. 

For each rotation rate Qq, a set of 2 000 images is ac- 
quired, covering at least 80 rotation periods. The sam- 
pling rate is synchronized with the platform rotation rate, 



with a number of images per rotation ranging from 24 
(for low Qq) to 9 (for large Qq). PIV fields are com- 
puted over successive images using 32 x 32 pixels inter- 
rogation windows with 50% overlap, leading to a spatial 
resolution of about 2 mm. This resolution is not enough 
to resolve the thickness of the Ekman boundary layers, 
6 ~ RE 1 ! 2 = 0.8 — 2.2 mm, but is appropriate for the 
large scales of the precession flow expected in the bulk. 

In view of the very low velocity expected for the pre- 
cession flow, the resolution of the velocity measurement 
is critical in our experiment. The characteristic velocities 
of the flow encountered in this work ranges from 0.01 to 
0.4 mm s _1 for Qq between 2 and 16 rpm. For the sam- 
pling rates considered here, these velocities correspond to 
a typical frame-by- frame particle displacement of 0.16 to 
2.6 pixels only. Although very weak, such displacement 
may actually be measured using PIV with sub-pixel inter- 
polation of the correlation peak. For interrogation win- 
dows of size 32 x 32 pixels, an accuracy of 0.05 pixel can 
be achieved using this technique [15 ,16 , yielding a signal- 
to- noise ratio ranging from 3 (low Qq) to 50 (large Qq). 

The orientation of the experiment with respect to the 
Earth rotation axis is monitored using a continuous laser 
beam aligned along the North- South direction and passing 
through the rotation axis of the sphere (see fig. [2]). The 
beam crosses the cubic glass tank and is therefore visible 
on the recorded images. The angle 6{t) between the South- 
North direction and the measurement fields (see fig. [2]) can 
be determined for each image with a precision better than 
±0.5°. 

Structure of the tilt-over flow. — We first show 
in fig. [3] the flow measured in the horizontal plane in the 
rotating frame for a rotation rate Qq = 6 rpm. This flow 
represents the departure between the total flow in the lab- 
oratory frame and the solid body rotation at Qq. In order 
to improve the quality of the velocity fields shown here, a 
phase average is performed over the velocity fields at the 
platform rotation rate Qq. This procedure allows to de- 
crease the broad-band PIV measurement noise by a factor 
iV 1 / 2 , where N is the number of recorded rotation periods 
(TV > 80). The spatial structure of the precession flow can 
finally be extracted with a signal-to-noise ratio of at least 
30 for all rotations rates. 

The four snapshots shown in fig. [3] are separated by a 
phase shift of 7r/4, with a phase origin chosen such that 
e x = e i (i-e. e y pointing to the North). In spite of the 
very weak velocity signal (of order of 0.04 mm s _1 , to be 
compared to the typical velocity of the sphere boundaries, 
fio R — 72 mm s _1 ), we clearly observe a well-defined 
flow pattern, which is rotating as a whole at the platform 
rotation rate but in the opposite direction. This weak flow 
is therefore stationary in the laboratory frame. Assuming 
that the total flow in the laboratory frame is a solid body 
rotation of vector cv slightly tilted with respect to Ho, the 
measured flow must be a solid body rotation of rotation 
vector uj' = uj — Qq. Since the measurement plane is 
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Fig. 3: Horizontal velocity fields measured in the rotating 
frame, in the off-centered horizontal plane at Zmes/R — 0.19 
for Q = 6 rpm [E — 1.2 x 10~ 4 ), with a phase shift of 7r/4 
between each image. The platform rotation is anticlockwise. 
The red arrows indicate the direction of the North at each 
time. Resolution of the velocity fields has been reduced by a 
factor 5 for better visibility. 

shifted at z mes /R ~ 0.19, the resulting horizontal velocity 
field must be uniform in the bulk, given by uj' x (z mes es), 
and rotating in the anticyclonic direction at frequency Qq , 
which is precisely what we observe. Snapshots at other 
values of Qo show essentially the same flow patterns. 

Measurements in the vertical plane, shown in fig.|4j con- 
firm this flow structure. In this configuration, the camera 
is now rotating around the vortex of quasi-horizontal ro- 
tation vector uj' stationary in the laboratory frame. The 
4 snapshots taken over half a rotation around the vortex 
actually show the following sequence: (a) anticlockwise, 
with uj 1 pointing towards the camera; (b) intermediate; 
(c) ascending, with u/ pointing to the left; (d) intermedi- 
ate. If the tilt-over flow were a pure solid-body rotation, 
the ascending flow in the snapshot (c) would be uniform, 
given by u/ x (y me s^y)^ which is approximately the case 
far from the boundaries. The wall region where the flow 
departs from a pure uniform flow has a thickness of order 
of 0.3 i?, which is much larger than the expected Ekman 
thickness E^ 2 R ~ 0.01 R. The tilt -over flow is therefore 
not exactly a pure solid body rotation, in agreement with 
numerical results obtained in a spherical shell with a very 
small stress- free inner solid core [T7]. Indeed, because of 
the breakdown of the Ekman layer at the so-called critical 
circles, a pure solid body rotation cannot be a uniformly 
valid solution [13]. 




Fig. 4: Vertical velocity fields measured in the rotating frame, 
in the off-centered vertical plane at y m es/R — 0.19 for Qo = 
6 rpm (E = 1.2 x 10 -4 ), with a phase shift of ty/4 between each 
image. The color maps the vertical velocity norm normalized 
by its maximum in each field. The phase origin is not the same 
as in fig. [3] 



Viscous prediction for the tilt-over flow forced by 
precession. — We compute here the rotation vector uj 
in the bulk of the fluid viewed from the precessing frame of 
reference (here the laboratory frame), following Refs. [T8l 
[T9] . The differential rotation between the fluid in the bulk 
rotating at cv and the sphere boundary rotating at Qq 
is matched across an Ekman boundary layer of typical 
thickness RE 1 ! 2 . We therefore assume E <C 1, such that 
a separation between a bulk flow and a thin boundary 
layer may be assumed. In the steady state, the viscous 
torque Y v exerted by the boundary layers on the fluid 
in the bulk is balanced by the Coriolis torque T c (note 
that the pressure torque is zero here because of spherical 
symmetry). This balance, projected along u> and along the 
two directions normal to u?, yields the following nonlinear 
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system of equations [5|fT9] . 

-^7= (c^3 cos A — co>2 sin A) = . . . 



(2) 



5 / 4 

x l/4 n 3/4 . x ^0 /o\ 



a 



^=uji cos A = A r r^ 4 c<j3^ 4 (^o — ^3)7 (4) 
V E 

where A r and A^ are respectively the non-dimensional vis- 
cous damping rate and viscous correction to the eigen- 
frequency of the tilt-over mode. Their values have been 
obtained by Greenspan [5] and completed by Zhang et 
at. [20], A r = —2.62 and A^ = 0.258. In presence of 
viscosity, the eigenfrequency ^0 of the inviscid tilt-over 
mode becomes VLq + (A$ + i\ r )\^^TtoTu [H[20], which 
means that, if the precession forcing is switched off, the 
tilt-over mode starts to rotate in the inertial frame at a 
frequency \i\fE\/Vto cJ, while exponentially decaying at a 
rate \X r \\^\^Tlocu. 

Equation © reflects the fact the work done per unit 
time by the viscous torque is zero, T v • uj = 0, since the 
work done by the Coriolis force is zero by definition. This 
equation, which can be recast into uj- (uj — fio) =0, simply 
expresses the so-called "no spin-up" condition, indicating 
that there is no differential rotation between the fluid and 
the sphere in the direction of the fluid rotation. This right 
angle between ui and uj' = uj — Qo indicates that the 
rotation rate \u>\ of the fluid is lower than Oq- 

If we further assume that the Poincare number CI p /CIq 
is small compared to E 1 ! 2 , the rotation vector u> is almost 
aligned with fio, and the system of equation ([2])-(|4j) can be 
simplified. More precisely, this regime applies for rotation 
rates 1]q > ^ c, with 



QpRs'mA 
\ r \fv 



(5) 



This condition is comfortably satisfied in the present ex- 
periments, with Qq,c — 5.2 x 10 -5 rpm. In this limit, the 
components of the tilt-over flow can be explicitly derived, 

1/2 

(6) 



0J2 



UJ 3 



Qp cos A 



A r 



CIqR 2 



A r 
fio- 



OJi, 



(7) 

(8) 



The horizontal projection of lo in the laboratory frame, 
u)h = oj\G\ + ^2^2, has therefore an amplitude 

1/2 / \2\ 1/2 

"I (1 + ^) > (9) 



I A, I 



which has indeed the expected form (pQ). Note that, in the 
limit considered here (Q p /Qo <^ V^), the horizontal pro- 
jection u?h measured in the experiment almost coincides 
with uj' . 
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Fig. 5: (Color online) Angle ip between the rotation vector of 
the tilt-over flow and ei (i.e. East) as a function of the rotation 
rate Qo in polar coordinates. Measurements are obtained in 
the horizontal (o) and the vertical (o) plane respectively. The 
continuous line shows the theoretical prediction cp = 174.35° 

Cod. 



In this limit, the angle ip between u>h and ei (the East 
direction) is constant, and given by 



(p = arctan 



— J = arctan ^ = 174.35°, (10) 



showing that u)h points almost to the West (along — ei), 
with a slight component to the North. Remarkably, this 
asymptotic angle obtained in the limit of large Qq is al- 
most perpendicular to the inviscid prediction of Poincare, 
for which Uh points to the North (i.e. cp = 90°). This 
indicates that, even for very low viscosity, the boundary 
layers have a critical influence on the tilt-over flow, pro- 
vided that Qp/Qo < E 1/2 . 

Comparison with the experimental tilt-over flow. 

— The rotation rate uoh of the tilt-over flow and its an- 
gle ip with the East have been systematically determined 
for f^o ranging from 2 to 16 rpm. These data have been 
extracted independently from the raw velocity fields mea- 
sured in the vertical and horizontal planes, and are com- 
pared here with the theoretical predictions ([9])- (fTU|) in 
figs. [5] and [6l 

Measurements of the vortex angle cp from the PIV data 
in the vertical plane have been obtained as follows: the 
horizontal vorticity, spatially averaged over a central re- 
gion of 50 mm radius, shows a harmonic oscillation at 
frequency Qq. At each period, the delay between the time 
t max of maximum vorticity (when u)h points to the cam- 
era) and the time at which the North-South laser beam is 
aligned with the camera axis is computed. Knowing the 
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Fig. 6: Horizontal rotation rate ujh of the tilt-over flow as a 
function of the rotation rate of the platform Qo, measured in 
the horizontal (o) and vertical (o) planes. The continuous line 
shows the prediction J9]). 

instantaneous angle 6(t) between the camera incidence and 
the South-North direction, we can simply deduce the an- 
gle of the vortex as (p = 0(t max ) + 90°. An independent 
estimate for ip has been determined from the data in the 
horizontal plane, by computing the time averaged (and 
spatially averaged over the region |r| < 50 mm) angle of 
the velocity with respect to the East direction ei. 

The rotation rate uoh of the horizontal component of the 
tilt-over flow has been determined from the vertical cuts 
as half the spatially averaged (over a central disk of radius 
50 mm) vorticity, measured at the times t max of maximum 
vorticity. uoh has also been determined independently from 
the horizontal cuts, as uoh = (|u^|)/^ mes , where (•) is an 
average over time and over the region |r| < 50 mm, and 
^mes is the height of the measurement plane. 

For both measurements in the horizontal and vertical 
planes, one value of ip and uoh is obtained at each rotation 
period. From this set, the average and standard devia- 
tion are computed over the 80 periods recorded for each 
rotation rate. In addition to the temporal fluctuations, 
the errorbars in figs. [5] and [6] also include the variations 
of (p and ujh when varying the radius of the averaging re- 
gion between 25 and 75 mm. For both quantities, the 
estimates determined from the two measurement planes 
closely agree, although data from the horizontal plane sys- 
tematically show a larger scatter. 

The vortex angle measured from both vertical and hor- 
izontal planes, cp ~ 173 ± 4° and 175 ± 11° respectively 
(fig. [5]), are in good agreement with the theoretical pre- 
diction ([TO]) 0. Similarly, the rotation rate ujh measured in 
both planes closely follow the prediction © to within 20% 
over the range Qq = 2 — 16 rpm (fig. [6]). The agreement 

1 A possible residual ellipticity of the sphere would lead to slightly 
different angles cp. Considering a prolate or an oblate spheroid, of 
ellipticity given by the maximum deviation of the radius of the sphere 
(R = 115 ±0.25 mm), yields predictions for ip between 170 and 180° 
for the range of Qo considered here, which is compatible with the 
present data. 



of uoh and ip with the theoretical predictions is remarkable 
in view of the very weak velocity signal, providing strong 
evidence that the weak secondary flow that we observe 
originates from the precession of the experiment by the 
Earth rotation. The magnitude of the secondary rotation 
lies in the range (1.5 — 3) x 10 _3 ^o 7 confirming that the 
rotation vector cv of the fluid is almost aligned with Ho, 
with a very weak angular departure of Uh/^o < 0.2°. 

Conclusion. — Measuring the influence of the Earth 
rotation at the laboratory scale is a technical challenge. 
In the fluid analogue of the Foucault pendulum presented 
in this letter, the very weak precession driven flow would 
have been impossible to detect directly from the labora- 
tory frame. Probing the flow in the rotating frame nat- 
urally subtracts the first-order rotation and allows us to 
detect this slight correction. We note that such residual 
tilt-over flow forced by the Earth rotation defines an irre- 
ducible background flow which should be present in every 
rotating fluid experiments, routinely used as models for 
geophysical and astrophysical flows in the laboratory. 
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